Derivation of the Gasser-Leutwyler Lagrangian from QCD 
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The normal part of the Gasser-Leutwyler formulation of the chiral Lagrangian is formally derived 
from the first principles of QCD without taking approximations. All the coefficients are expressed 
in terms of certain Green's functions in QCD, which can be regarded as the fundamental QCD 
definitions of the normal part of the coefficients. The approximate values of the claculated coefficients 
are also presented. 
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Chiral Lagrangian (CL) ^ Q, especially the three- 
flavor formulation given by Gasser and Leutwyler H, 
has been widely used in the study of low energy hadron 
physics. Studies on understanding the relation between 
the CL and the underlying theory of QCD will be very 
helpful for making the theory more predictive. There 
arepapers studying approximate formulae for the CL 
[^-^ either based on certain assumptions or considering 
the anomaly contributions with certain approximations 
from the beginning. Further improvements are certainly 
needed. Actually, the study can be divided into two 
steps: (i) formally deriving the CL from QCD giving the 
fundamental QCD definitions of the coefficients, (ii) cal- 
culating the the values of the coefficients from the QCD 
definitions with certain approximations. This paper is 
mainly devoted to step (i) in which we derive the fun- 
damental QCD formulae from the normal part contribu- 
tions of the theory without taking approximations, which 
is different from and complimentary to the anomaly part 
in Ref. pj. We shall also present the values of our ap- 
proximately calculated coefficients [step (ii)]. 

Let A"-^ be the gluon field, i/'S'' and *J^'' be, respec- 
tively, the light and heavy quark fields with color index 
a, Lorentz spinor index 77, light flavor index a and heavy 
flavor index a. Following Gasser and Leutwyler, we start 
from the generating functional M 



ator ■0'"» (a;)-0'"'(a;) and integrate out all the remaining de- 
grees of freedom. For this purpose, we consider the scalar 
and pseudoscalar degrees of freedom of {p'"' (a;)i/'"'' (x) and 
make the following decomposition 



Z[J]^ 



where £ is the QCD Lagrangian with the gauge-fixing 
term and the Fadeev-Popov determinant. The local ex- 
ternal sources Ja-p can be decomposed as 

J{x) = -s{x) + ipix)-f5 +'^ (x) + (^ (2^)75, (2) 

in which the light quark masses have been absorbed into 
s. Since the contributions from the anomaly to the CL 
has been studied in Ref. g , our aim is to study the com- 
plete normal part contributions, and thus we simply take 
the 0- vacuum parameter 6 — 0. 

What we need to do is to consistently extract the pseu- 
doscalar degree of freedom from the local composite oper- 



^'H^){i5)c,r"ix) 






fi'Vfi't)"''(a;) 
-n"'an"')''\x) 



(3) 



where the modular degree of freedom is described by 
an hermitian field (j{x) [<7^(x) — cr{x)], and the phase 
degree of freedom is described by a unitary field fl' 
[Q'^{x)n'{x) — 1]. As usual, we can define U'{x) = 
n''^{x) which contains a U{1) factor such that detU'{x) = 
e"' ''^'' {^'(x) is real). We can extract out the [/(I) fac- 
tor and define a field U{x) as U'{x) = e~ U{x) with 
c\ctU{x) = 1. Then we can define a new field fl and 
decompose U into 



u{x) = n^{x), 



(4) 



as in the literature. This U{x), as the desired repre- 
sentation of SU{Nf)fi X SU{Nf)L, will be the nonlinear 
realization of the meson field in the CL. It is straightfor- 
ward to eliminate a from the two equations in (H) and 
get the following relations. 



e ■'"/ nhTi[PRi;^]n'< = e"f ntTi[PL^ip]n, 



PV2?V;P*I?i'PA^e'/^"=^^''^•'^'*•*"^^)+'^'''^^ (1) 



o2iiJ' 



det[tri[PR^^]]/det[tvi[PLij^]], 



(5) 



(6) 



in which all the fields are at the same space-time point 
X, and tr/ is the trace for the spinor index. Eqs.(p[)- 
(||), especially (||), describe our idea of localization in the 
operator formalism. 

To realize this idea in the functional integration, we 
need a technique to integrate in this information to (|l|). 
For this purpose, we start from the following functional 
identity for an operator O satisfying dot O = dot O^ |g| , 

fvu 5{u^u - i)5{detu - 1) j'p] 5{no^n - n^on'') 

= const , (7) 



in which f VU d{U^U — l)5{detU — 1) is an invariant 
integration measure and the function J-'[0] is defined as 



T[0] 
With 



i- = detO f Va 6{OlO - a^a)6{a - a^) 



■ i9'(x) 



0{x)=e "/ traPfl^(x)^(x)], 



(8) 



(9) 



(R) serves as the functional expression reflecting (H). 
Inserting ^ and (g) into (|l), we obtain 

Z[J] = / Vi^V^jV^V^VAf^VU 5{U^U - \)5{deiU - 1) 



X(5 e'"/mr/[V'LVi?J^ 



'«f ot 



nhri[^R^L]n^ 



X exp {«r/[^] + i I d''x{C{4>, i>, *, §, A^) + t/; JV-}} 



DC/ (5(C/^C/ - l)5{deiU - 1) e*^'=^^['^'"'l, 



(10) 



where 



^ ^[0(x)] 
= n { [det{^tr,e[V'i?ViL]}det{^trz,[^LV^ii]}] ^ / Da 

a; "^ ^ 

xS{ — tTic{^Ri,L)tM^L^R) - CrV)(5(CT - (7^)}, (11) 

and we have introduced the effective action 



N, 



+i / d''x{/:(7/., i/;, *, *, A^,) + ipjilj} 



(12) 



to formaUy express the integration over the quark and 
gluon fields. 

To work out the C/(a;)-dependence of Sgff, we make 
the following chiral rotation 

Mx) = [n{x)PR + n\x)PL\[J{x) + ip] 
yi[n{x)PR + n\x)PLl 

ijn(x) = [n\x)PR + n{x)PL\ilAx), 

Mx) = ^p{x)[n\x)PR + nix)PL]. (13) 

Here we have denoted the rotated quantities by a sub- 
script H. The fields C/, >1^,^ and ^ are unchanged 
under the rotation. It is easiy to see from (H) that 
i?Q ~ "d' ■ In terms of the rotated fields, the (5- function 



in (12) is (5 e "ftnlML^i 



'n,R\ 



''fiul'ipn.R'ipn.^ 



so that the explicit [/-dependence of the theory is from 
the source term and the ^-function in (p^). The re- 
maining part in (03) is invariant under rotation (n3) 
since the Lagrangian £, together with the source-term, 
is chirally invariant, and from (O) we can see that 
p [^] — r [ '^y " ] . Since the integration measure VipV^ 
is not invariant under the chiral rotation, it will cause cer- 
tain anomaly terms in S'e// i-e. / Dt/'Dt/'D^D'I'DA^ = 
J VtPnV^^V^V^VAf, el^'^^o^aly terms] ^ ^pj^^g^ making 
a chiral rotation of -ip and ■0, we have 

+i / d'^x[C{ip, -0, ^, ^, A^) + i/J JnV'] + anomaly terms} 



ViJV^|JV^V■i'VAf, (5( !/;"( 



ijj'tjj 






- 75 cos 



Nf 



)i/ 



exp {ir/[^] + ij d^x[C{ip, 0, *, *, A^) + i^ifn 

■d' 
+(^n75 - sn ~ pn tan -rp)V''] + anomaly terms} . (14) 



In (14) and all the later formulae, we simply use the sym- 



bol tp (tp) as the short notations oitpQ {tpn) for shortening 
the formulae. This does not make confusion since they 
are integration variables). We see that, in (M), the U- 
dependence of Seff comes only from the rotated sources. 
So that the rotation simplifies the [/-dependence of Seff- 
Eq.([l4|) shows that S'e//[1, Jq] is the QCD generating 
functional for the rotated sources with the — zt/)"-;/)'' sin -^ 
+'tp°''y5ip^ cos 4- degree of freedom frozen. After making 



Nf 



i3' 



a further [/^(l) rotation of the ip and ip, ^ can be ro- 
tated away and the frozen degree of freedom just becomes 
the pesudoscalar degree of freedom 'tjj°'^5ijj'' as it should 
be since this degree of freedom is already included in the 
[/-field. The automatic occurance of this frozen degree 
of freedom implies that our way of extracting the [/-field 
degree of freedom is really consistent. 

We first consider the p^-order terms. Expanding ([l4| ) 
up to the order of p^ , and taking account of translational 
invariance, parity and flavor conservations, we obtain 



Seff[l:Jn] 
^ fdxtr 



F^ / dx tT[al + Bosn] 



(15) 



where V^ is the covariant derivative related to the ex- 



ternal sources deflned in Ref. ||], x = 2i?o(s -I- ip), and 
the coefficients Fq and Bq are c 
following QCD Green's functions as 



the coefficients Fq and Bq are defined in terms of the 



fo^Bo^-^(^V'), 



(16) 



and 



_ /Cii ^8 r -Sq/Ci + NfK-ia + BqICis 



F, 



° ~ 8(iV2 - 1) 



dx 



i, io)ri5r{o)rixh^.i5nx) 



16^2' 



i8 = 



1 /r-va 



TVf 



--^6/ 



[^ (0)7''75V'"(0)][V (xh^l^rix)] 



(17) 



C i^8 



il 



WBo 



15 



16^2 



IC2+IC 



13 



, i^l 



.L9 = 
A^2 — /Cl3 



4/C13 - ^14 



In (nq) and (O), the vacuum expectation values are de- tt _ Bf^K-i + JC7 + ICg -Bq/Cis 
fined as, for an operator O, " "^ 



8Bi 



(20) 



0)^I^^O 



The coefBcients /Ci, • • • /C15 are defined by the following 
r pu ' integrations of the QCD Green's functions 



xe 



i^il^]+i f <lxC{jP,i^,'i/,^,Ai_,) 



,/Ci 



c 

d j:hc 



(18) = [(-1 - /C2)(g^^'.9''''' + g^'^'g"'"') + 21C2g^"'g''''''W6' 



and {• ■ ■)c denotes the connected part of (• ■ •). The in- 
tegrand in (lla) is just the p2_order terms in the Gasser- 
Leutwyler Lagrangian ||]. 

The p''-order terms can be worked out along the same 
line. With the help of the p2_order equation of motion, we x[i/)'^(y)7'^75-(/;^(2/)][-(/'^(2;)7'^75i/)''(z)] 
obtain the normal part contributions (ignoring anomaly 
contributions) to the p^-order terms, 



2^ ^ dxdydz(\ris>)ri5i'\mr{x)Yi5^\x)] 



Seff[l,JQ\ 



0{p* 



dx tr 



6 



S'=fS''^S3''[-{g^'•'g^'' + <?^^g''^)/C3 + g'^^g'^'^ICi] 



/Ci [d^a: 



Ml2 



+jchjg/^e6r ^^M^^Ak _^ g^^>^ g^'')]C:i + g'^'g^^ICi 



-/C2(d^a^ - d" a'^-{){dfj,an^i, - d^an,^) + /Caion]^ 
-|-/C4a[^aQaa,^ao,,. + /C5aQtr[a^] + /CeanOnti'I'^si,^^^^,!.] 
+/C7S^ + /C8Sntr[sn] + /Cgp^ + /CioPotr[pn] + /CnSna^ 



-K.ibPii.d^ a,n„ 



+5- 



■o.f 



S<^dgchggb,f^g^Xg^n ^ gt^'^g^^)jc^ + g^'^g^'^ICi 



^gehggdgcb^f^g^Xg^^ + gt"'g^'^)]C3 + g'^'g^'^ICi] 



-6' 



ah 



53dScfgeb^_(^g„^.gy\ _^ gt^>^g^^)]C^ + .g''''g^'=/C4] 



dx 



Li[tr(V^C/tV^C/)]2 






+L2tr[V^C/^V^C/]tr[V^C/tV'';7] 

+L3tr[(V''t/tV^C/)2] + L4tr[V''C/tV^C/]tr[x^C/ + U^] +5"^r^<5=''5«'^[5^^g'"'2/C5 + (g^''^^'^ + 5^''5''')^6] 

+L5tr[V''t/tV^C/(x^C/ + U\)] + i6[tr(x^C/ + C/^x)]^ 



+L7[tr(x^C/ - t/^x)]' + i8tr[x^ C/x^C/ + xU^ xU^] 



+5<'h^gb^cfged^gt.KgU\2^^ + (g^^g^'^ + 5^^g'"')/C6] ^ + • • 



- / dx{ [r{o)^'m[r{x)i^''{x)] 



ab xcd 



+i?2tr[x^x] 



(19) 



where d^aj^ ^ d^a^^ - tv^a^^ + za^«^, V^" ^ df^v^, 
d^v^ - »<«o + ivnVn^ and 



'- J dx(^[r{0)^\0)] tan^[V;^(a;)V/(x)] tan '^'^''^ 



TV 



/ /c 



= ICa5'"^S^'' + ICwS''''S'"^ 



Fi 



/C4 + 2/C5 + 2/C13 — ^14 



1 



^2 = 



32 

/C4 + /Ce + 2/C13 — A^i4 



dxdyn^''(0)V'(0)][V;^(x)7'^75^''(x)] 



x['0''(2/)7m75'0 (2/)] 



16 



_ /C3 - 2/C4 - 6/C13 + 3/C14 _ /C12 

-f^3 — 77 : J^4 — 



llCn{S'"^S''fS'''' + 5''f6'"^5''^) + /Ci2<5"''<5"^<5^'^ + • • • , 



16 



I6F0 



/c 



13 



576(7V? - 1) 



dx 



(,&5m'^5'm''^' ^.9/j/j'5i'f' J'*' -'' 



x[ [^ {o)rrm [rixh'-rix)] 



Nf 



[r{o)rrm[i^ i^h^rix)] 



ICi 



ICl5 






c\ rriU'uXK 

^9tJ.iy9\K.^A 



gtj.\gvt,±A 



o rrifiiyXK r% miivXK, . mUyXK. 



(Tx xf" 



V'" (0)1/^^(0) tan 



7?'(0) 



4(7V; - 1) 

xV;^a;)7,75^°(a^)) - ^^^^'^(0)^(0) tan ^ 



xV'''(a;)7M75V'''(a;) 



with Ta, Tg defined as 



(21) 



-- / dxdyx-{ [r{0h''ri0)Mixh''l5rix)] 



xifivh^isi^Hy)] 



^adccf cebrj^l-ii^Xn , :af ged ^chrpiivXK. 



(221 



in which the unwritten terms are those which do not 
contribute in ([l9|). The integrand in ( [l9| ) is just the p^- 
order terms in the Gasser-Leutwyler Lagrangian. 

Eqs.([l6l), (I?!), and (|2|)-(|2|) can be regarded as the 
fundamental QCD definitions of the normal parts of 
Fq,Bq,Li,---, Lio, Hi and H2. In principle, the related 
Green's functions can be calculated from lattice QCD. 

As our first calculated result, we present here the val- 
ues of the coefficients calculated from solving the equa- 
tions for the related Green's functions in a crude ap- 
proximation. We consider only the contributions from 
the two-point Green's functions, and take the approxi- 
mations of the large- TVc limit and the leading order in 
dynamical perturbation theory Hj4|. In these approx- 



imations, our formula ( [l7| ) for F^ reduces to the well 
known Pagels-Stokar formula m. In the present ap- 
proximation, the QCD effective coupling constant is ap- 
proximately as{p'^) ~ 5s/(47r[l + Hcip'^)]) |^, where 
IIg is the gluon self-energy. In the evaluation of the 
Green's functions, we further take the approximation 
<^s[{p — q)'^] ~ as[max{p'^ , q"^)] to reduce the integral 
equations to differential equations. In the calculation, we 
take the modified one- loop formula for aj, i.e. Usip^) — 

(33-N 'iinf^VA'^ +tT B ' ^^'^'^ rcnormalization schemes 

cannot be distinguished in the one-loop formula and we 
take Nf = 3, the value of Aqcd is not to be compared 
with the experimental value of A— — . 

^ MS 



In the numerical calculations, we take Fq = 87 MeV 
[ [lO[ as input by which r is related to Aqcd- We present 
here, in Table 1, the numerical results of the complete 
coefficients Li, ■ ■ ■ , Liq, Hi and H2 (the normal part plus 
the anomaly contribution 0) for Aqcd = 600, 700, 800 
MeV (see Ref. |ll for details). 

Table 1. Values of the complete p''-order coefficients (in 

10"'') for various values of Aqcd, and the comparison with 

the values determined by experiments from Ref. pi. 



Complete 
Coefficients 


Aqcd (MeV) 


Expt. 


600 


700 


800 


Li 


1.0 


1.0 


1.0 


0.9± 0.3 


L2 


2.1 


2.1 


2.1 


1.7± 0.7 


L3 


-3.1 


-2.5 


-2.0 


-4.4± 2.5 


U 











0± 0.5 


L, 


4.0 


4.7 


5.4 


2.2± 0.5 


Le 











0±0.3 


L7 


-0.27 


-0.41 


-0.50 


-0.4± 0.15 


Ls 


0.006 


0.44 


0.69 


1.1±0.3 


L9 


6.3 


6.3 


6.3 


7.4± 0.7 


Lw 


-2.4 


-2.1 


-2.1 


-6.0± 0.7 


Hi 


1.2 


1.1 


1.0 




H2 


-0.01 


-0.89 


-1.4 





We see from Table 1 that all results are of the right sign 
and most of them are very close to the experimental val- 
ues except L5, Lg and Liq. Existing certain deviation 
is understandable since the present approximations are 
rather crude. Improvement of the approach is in progress. 
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